Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 30 July 2012 (MN M£K style file v2.2) 



The effect of primordial black holes on 21 cm fluctuations 

Hiroyuki Tashiro 1 and Naoshi Sugiyama 2 ' 3 

1 Physics Department, Arizona State University, Tempe, AZ 85287, USA 

2 Department of Physics and Astrophysics & Kobayashi-Maskawa Institute for the Origin of Particles and the Universe, 
Nagoya University, Chikusa, Nagoya 464-8602, Japan 

3 Institute for the Physics and Mathematics of the Universe (IPMU), The University of Tokyo, Kashiwa, Chiba, 277-8568, Japan 



(N 

o 



30 July 2012 



ABSTRACT 



(N 

o 

u 

Of 

6 



> 
in 
o 

\o 

\>- 
o 

(N 



The 21 cm signal produced by non-evaporating primordial black holes (PBHs) 
is investigated. X-ray photons emitted by the accretion to a PBH ionize and heat 
intcrgalactic medium (IGM) gas near a PBH. Using a simple analytic model, we 
show that this x-ray heating can produce an observable differential 21 cm brightness 
temperature. The region of the observable 21 cm brightness temperature can extend 
to 1-10 Mpc comoving distance from a PBH with depending on the PBH mass. The 
angular power spectrum of 21 cm fluctuations due to PBHs is also calculated. The 
peak position of the angular spectrum depends on the PBH mass, while the amplitude 
is independent. Comparing it with the angular power spectrum caused by primordial 
density fluctuations, it is found that both of them become comparable if ripBH = 
l(r 12 (A//10 3 Af Q )-°- 2 at z = 30 and 10- 13 (Af/10 3 Af )-°- 2 at z = 20 for the PBH 
mass from 10 M to 10 8 M©. Finally we find that the Square Kilometer Array can 
detect the signal due to PBHs up to fi PB H = 10 _9 (M/10 3 Af©)" 2 at z = 30 and 



10- 10 (A//10 3 MqY 



at z = 20 for PBHs with mass from 10 2 Mr, to 10 8 Mr, 
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1 INTRODUCTION 



It is believed that primordial black holes (PBHs) formed in the early Universe I Carr fc Hawkinsl 1974 : Carr 19751 ). Although 



there is no direct evidence of the existence of PBHs, PBHs are attracting attention as the way of constraining physics in the 
early Universe. In particular, one of the main generation mechanisms of PBHs is a gravitational collapse of an overdensity 
region at the horizon scale when the amplitude of overdensity exceeded a critical threshold. Therefore the resultant mass 
functio n and the abunda nce of PBHs depend on the amplitude of primordial density fluctuations at the horizon crossing 
epoch ( Green et al. 20041 ). The PBH abundance is expected to be a probe of primordial density fluctuations on small scales 



which cannot be accessible by CMB or LSS observations. 

Constraints on the abundance of PBHs have been studied in many papers and continue to be updated ( Carr et al. 2010l ). 
PBHs with mass smaller than 10 15 g have evaporated by the present epoch because the evaporation time scale by Hawking 
radiation is less than the Hubble ti me scale today jHawkinj 19741 ). However evaporation of PBHs g enerates additiona l 

'A 



entropy in the Universe after inflation llze l'dovich fc Star obinskiil 19 76|). affects on b ig bang nucleosynthesis ( Vainer et al.lll97^ : 



Vainer fc Naselskiill 19781: iMivama fc Satolll978l : IZeldovich et al.lll977l : lLindlevlll980l ) and distorts the CMB bl ackbody spectrum 
lTas!in^^r*sTgivEUTia 20081). PB H evaporation may also produce the observable gamma-ray background l|Page fc Hawkinsl 



19761 ; iMacGibbon fc Car? 199 ll ). According to measurements of these cosmological phenomena, there are strong constraints 



on the abundance of PBHs with mass smaller than 10 15 g. 

PBHs with mass larger than 10 15 g survive in the present Universe. One of the constraints on such PBHs can be set 
from the fact that the current density parameter of PBHs, JIpbh, cannot exceed the cold dark matter density parameter 
observed at the present epoch, fic- Conventionally, the constraint on P BH abunda nce is given by P(M) which is the fraction 
of the regions of mass M collapsing into PBHs at the formation epoch l|Carrlll975l ). The constraint on the de nsity parameter 



of PB Hs today, li PBH < ^c, implies j3 < 2 x 10~ ls (M/10 15 g) 1/i from WMAP 7-year dat a, i.e., f2 c = 0-22 ( Komats u et al 



20 111 ). Microlensing observations also constrain the abundance of non-evaporating PBHs IjAlcock et al.ll200ll ). iRicotti et al 
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(2008) have obtained the constraint on PBHs with mass larger than 0.1 M@, investigating the effects of such PBHs on cosmic 
reionization and CMB temperature anisotropies and distortions. The future gravitational w ave observations are also expected 
to provide a probe of the massive PBH abundance (jloka et al.lll999l ; llnoue fc Tanakall2003T ). 

In this paper, we evaluate the 21 cm brightness temperat ure produced by PBH s and study the potential of 21 cm 
observations to give the constraint on the abundance of PBHs. Mack fc Wesley 1 20081 ) have investigated the signatures of 
evaporating PBHs in 21 cm brightness temperature. Accordingly, they have concentrated on PBHs whose mass range is 
5 x 10 13 g < A4"pbh < 10 17 g. On the contrary, here, we focus on non-evaporating PBHs with mass much larger than 10 15 g. 
It has been shown that PBHs with large mass could produce X-ray and UV photons through the accretion of matter to 



-fr 

PBHs and these photons heat up and ionize intergalactic medium (IGM) (Carr 1981; IGnedin et alj 1995; Mi ller fc Ostriker 
2001 ; Ricotti et al. 2008). The heated and ionized IGM gas may produce the observable deviation of the 21 cm brightness 



temperature from the background. Our aim in this paper is to evaluate this deviation and to discuss the potential of the 
21 cm observation to give the constraint on the non-evaporating PBH abundance. 

The paper is organized as follows. In Section 2, using a simple model of X-ray photon flux due to the accretion to a 
PBH, we evaluate the ionization and temperature profile near a PBH. In section 3, we calculate the spin temperature and the 
brightness temperature induced by a PBH. In section 4, the angular power spectrum of 21 cm fluctuations due to PBHs are 
evaluated. Section 5 is devoted to the conclusion. Throughout the pape Throughout this paper, we use parameters for a flat 
ACDM model: h = 0.7 (Hp = h x 100 km/s/Mpc), S1 B = 0.05 and Q M = 0.26. These parameters are consistent with WMAP 
results IjKomatsu et al.ll201ll ). 



2 IONIZATION AND HEATING OF IGM BY A PBH 



The accretion disk of a non-evaporating PBH, whose mass is larger than 10 15 g, can be a source of X-ray photons. For example, 
the non-evaporating PBHs are one of the candidates for super massive black hole seeds. PBHs whose masses exceed 10 5 Mq 
cannot be directly formed through gravitational collapse since the time scale of the collapse becomes longer than cosmological 
time. Accordingly on ly the accretion after the formation, which induces X-ray photon emission, makes possible to form such 
massive black holes (jPuchtingl bool iKhlopov et al.ll2005l) . It is difficult to theoretically predict the accurate X-ray photon 
spectrum from PBHs, because the X-ray spectrum depends on the detailed condition of the accretion and the environment of 
PBHs such as the amount of neutral hydrogens. Hence, for si mplicity, we assume that the PBH accretion powers a miniquasar 
with the power-law spectrum of X-ray photons, according to lKuhlen fc Madaul l|2005l ). 

F(E)^AE' 1 s~\ (1) 

where A is set as the tenth of the Eddington luminosity. Here we consider the range of the photon energy E from 200 eV 
to 100 keV since we assume that emit ted photons whose e nergy are lower than 200 eV are immediately absorbed by the 
surrounding gas of the PBH. Following IZaroubi et alJ l|2007h . we evaluate the ionization and heating of IGM due to massive 
PBHs in this section. 

The number density of photons per unit time per unit area at distance r from the source is given by 

A 



N{E-r) = e~ T{E '' r) 



T (E;r) = 



(Anr 2 ) ' 
nHXH&(E)dr, 



(2) 



(3) 



where xh is the hydrogen neutral fraction, uh is the mean number density of hydrogen at a redshift z and cr(E) is the 
absorption cross-section per hydrogen atom. In order to take into account t he contribution from helium atoms as well as 
hydrogen atoms, we adopt the fitting formula bv lZdziarski fc Svensson jl989l ). 



a = 4.25 x 10" 



E 



250 eV 



2.65 for E < 250 eV, 
3.30 for E > 250 eV. 



Using the function N(E\r), we can write the ionization rate per hydrogen atom at the distance r from the source as 



r(r)= / a(E)M(E;r) 
Je 



E A dE 



(4) 



(5) 



where the term, (E/Eo)<f)(E, x e ), is introduced to consider the secondar y ionization by the p hotoelectrons produced by 
energetic photons ( E > 100 eV). We app l y the fitting formula of 4> by iDiikstra et al.l (2004) for the low energy region 
E < 0.5 keV and bv lShull fc van Steenberd (|l985l ) for the high energy region E > 0.5 keV. 

The neutral fraction of hydrogen is obtained by solving the equation of the ionization-recombination equilibrium 



OLHn H (l — xh) 
where a.n is the recombination cross-section; ani 



T(r) uhxh 



(6) 
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comoving radius [Mpc] comoving radius [Mpc] 

Figure 1. The neutral fraction of hydrogen as a function of comoving distance from the source. The left panel shows the dependence 
on mass at z = 30. The lines represents the neutral fraction with masses ranging from 10 Mq to 10 7 Mq from left to right. The right 
panel shows the dependence on redshift for a PBH with M = 10 3 Mq. The dottcd-dashed, dotted, solid and dashed lines represent the 
neutral fraction at z = 50, z = 40, z = 30 and z = 20, respectively. 



Solving Eq. we show the neutral fraction of hydrogen for different masses and for different redshifts in the left and 
right panels of Fig.[T] respectively. Here we assume that the density around a PBH is same as the mean density of the Universe. 
The comoving radius of the ionization sphere is made large by massive PBHs. Increasing mass of PBHs means enhancing the 
ionization photons, because we assume that the flux of the ionization photons is proportional to a tenth of the Eddington 
luminosity. As the Universe evolves, the comoving radius of the ionization sphere increases slowly. 

Next, we evaluate the kinetic temperature of IGM around a PBH. The heating rate per unit volume per unit time at the 
distance r form the source is obtained by considering the photons absorbed by the IGM at r, 

U(r) = fn H x H (r) a{E)M{E;r)dE (7) 
Je 

where / is the fraction of the absorbed photon energy through the collisional excitations of the IGM. Shull fc van Steenber d 
! 1985t ) provided a simple fitting formula / = C \l - (1 - x a ) b ^ , where C = 0.9771, a = 0.2663, b = 1.3163 andi^l-in is 
the ionized fraction. 

The kinetic temperature of the IGM at the distance r, Tk(r), is determined by the balance between the heating and the 
Compton cooling due to CMB photons. 

W(r) = f^r 7 4 (l - x„){T k {r) - T 7 ) + 2HT k {r), (8) 

where or is the cross section for the Compton scattering and T 7 is CMB temperature. Here we also take into account the 
cooling by the expansion of the Universe. 

Fig. [5] shows the gas temperature profiles for different masses in the left panel and for different redshifts in the right 
panel. Here we add the background temperature to the kinetic temperature in order to match both temperatures at a large 
distance from a PBH. Near the source, the temperature is determined by the heating rate and the Compton cooling rate. 
With increasing the distance from the source, neutral fraction grows and the optical depth r becomes larger. Accordingly, 
the number density of photons damps as shown in Eq. J2}. As a result, the temperature starts to decrease rapidly. Because 
the Compton cooling depends on the number of free electrons, this cooling becomes ineffective at the distance where the 
neutral fraction of hydrogen becomes almost unity. For example, this scale corresponds to 0.1 comoving Mpc for a PBH with 
M = 10 3 Mq at z = 30. Beyond this point, the temperature mildly decreases due to the cooling of the cosmic expansion. The 
kinetic temperature at the inner side is independent on a PBH mass. However the region of the high temperature becomes 
larger as the PBH mass increases. In the right panel of Fig. [2] the redshift dependence of the kinetic temperature is shown. 
The larger the neutral hydrogen density is, the larger the heating efficiency becomes as in Eq. 0. Therefore the heated 
temperature becomes high as the redshift increases. 



3 21 CM BRIGHTNESS TEMPERATURE DUE TO A PBH 

As shown in the previous section, a PBH ionizes and heats the surrounding region by the X-ray emissions. These ionization 
and heating leave an observable signature as the differential brightness temperature of the 21 cm intensity relative to the 
CMB temperature. 

The 21 cm intensity depends on the ratio between the number density of neutral hydrogen in the excitation and ground 
states in the hyperfine structure. This ratio can be quantified by the spin temperature T s . The spin temperature is determined 
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comoving radius [Mpc] comoving radius [Mpc] 

Figure 2. The kinetic temperature of hydrogen gas as the function of the comoving radius. The left panel shows the dependence on 
a PBH's mass at z = 30. The right panel shows the dependence on the redshift for a PBH with M = 10 3 Mq. In both panels, the 
representations of lines are same as in Fig. [T] 



by the balance between three processes: absorption o f CMB photons, collisional excitation in the kinetic temperature T k , 
and Lyman a pumping dWouthuvsenlll952 ; Field 19581 ). In the steady state approximation between these processes, the spin 



temperature is obtained by (jFieldll 19581 ) 



T» + Tj + y k T k + y a T k 



O) 



1 + yk + y a 

where T, is 0.068 K which is the temperature corresponding to the energy difference between the levels in the hyperfine 
structure, and yk and y a are the kinetic and Lyman-a coupling efficiencies, respectively. The kinetic efficiency is given by 

T, 



Vk 



■ (Ch + C e + C p ) 



(10) 



AloTkin 

where Aio is the Einstein spontaneous emission rate coefficient, Aiq = 2.9 x 10 -15 s _1 . The terms, Ch, C e and C p , represent 
the de-excita t ion ra tes due to neutral hydrogen, electrons and protons, respectively. Here we use the fitting formula by 
Kuhlen et all fcood ). 

CH=n H K, C e =n e 'y e , C p = 3.2n p K, (11) 

where n e and n p are the electron and proton number densities, respectively. In Eq. (| 1 1 f) . k is the effective single-atom rate 
coefficient, 



3.1 x 10~ 11 Tz ff r fc °' 357 exp(-32/Tfc) cmV 1 , 



(12) 

and 7e is given by log ( 7e / lcm 3 s" 1 ) = -9.607+0.5 log T k xexp(-(log T k ) 4S /1800) for T k < 10 4 K and 7e = j e {T k = 10 4 K) 
for T k > 10 4 K. 

The Lyman-a coupling efficiency is given by Field ( 19581 ) , 

_ 16n 2 T t e 2 f 12 J 
Va - 27A 10 T k m e c ' 

where /12 = 0.416 is the oscillator strength of the Lyman a transition, and e and m e are the electron charge and mass, 
respectively. In Eq. (|13[) . Jo is the flux of the L yman a photons due to collisional excitations. At the distance r from the 
source, Jo can be written as (jZaroubi et al.ll2007l) 

! C 



Mr) = 



-n H x H 



(r) / a(EW(E;r)^, 

J En nVa 



(14) 



where <j> a is the fraction of the absorbed energy going into the collisional excitation of Lyman a. lShull fc van Steenberd j 19851 ) 
gave the following analytical form, 

^^0.48(l-(l- a;ff ) a27 ) 1 ' 52 . (15) 

Now we can calculate the spin temperature in case existing a PBH according to the results in the previous section. We 
show the results in Fig. [3] The left panel shows the spin temperature for different masses of a PBH at z — 30 and the right 
panel represents the spin temperature for different redshifts for M = 10 3 Mq. 

In the highly ionized region, the profile of the spin temperature is flat and does not depend on PBH's mass as shown 
in the left panel of Fig. [3] Since it is highly ionized, the spin temperature is determined by the CMB temperature and the 
term of the kinetic efficiency. As the neutral fraction goes up with increasing the distance, however, the Lyman a efficiency 
y a becomes larger so that the term of Lyman a coupling begins to be effective. Peaks of T s in the small radius for small 
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Figure 3. The spin temperature as the function of the comoving radius. The left panel shows the dependence on PBH's mass at z = 30. 
The right panel shows the dependence on the redshift for a PBH with M = 10 3 Mq. In both panels, the representations of lines are 
same as in Fig.[T] 



PBH masses shown in the left panel are caused due to this effect. It should note that the peak values are also independent on 
mass since T 3 approaches to Tk for a large value of y a . Because the flux of Lyman-a quickly decreases as the radius increases, 
the spin temperature also drops. Finally, the spin temperature settles down to the one determined by the background CMB 
and kinetic temperatures. The larger the PBH's mass is, the larger the region where the spin temperature is higher than the 
background one becomes. 

Because the kinetic temperature strongly depends on the redshift as shown in the right panel of Fig.O the amplitude of the 
spin temperature also has the strong dependence on the redshift as shown in the right panel of Fig. [3] The spin temperature 
is almost same as the kinetic temperature in the ionized region at high redshifts (z > 40), since the kinetic efficiency yk 
dominates other contributions in Eq. ((9|. In low redshifts (z < 30), we find the region where the spin temperature is lower 
than the background one at the edge of the heated region. In these redshifts the background spin temperature is almost 
the CMB temperature, while the background spin temperature is between the CMB and kinetic temperatures in redshifts, 
30 < z < 50. Since the kinetic temperature is lower than the CMB temperature at the edge, the Lyman-a coupling is strong 
and draws the spin temperature toward the kinetic temperature. Theref ore, there exist s a reg ion where the spin temperature 
becomes lower than the background one. This tendency appears even in lZaroubi et al. (2007). 

The differential br ightness temperature from the CMB temperature for a given spin temperature T s is obtained by 
i Ciardi fc Madaull2003r i 

,lV2 



ST b = (20 mK) (1 



0.0223 



1 + z 
10 



0.24 \ 

Om / 



(16) 



where 5 is the density contrast. Here we assume 8 — as mentioned in section 2. In Fig. [4] STb is shown as a function of 
the comoving radial distance from a PBH. The left panel is for different masses of a PBH at z — 30 and the right panel is 
for different redshifts for M = 10 3 M Q . The ST b is almost zero near the source, because hydrogen in such regions are totally 
ionized. As the neutral fraction increases, the brightness temperature also grows. The peak amplitudes of 8Tb are independent 
on the PBH mass for M < 10 4 Mq. This is because the spin temperature is much larger than the CMB temperature, T B ^> T 7 . 
Accordingly, from Eq. (I16|l . STb becomes independent on T s and only depends on (l + z) 1//2 . This redshift dependence is shown 
in the right panel of Fig. 2] On the other hand, the peak amplitude becomes smaller for a mass of the PBH higher than 10 s M Q 
in our model. 

As in the case of T s , the region where the differential brightness temperature is below the background value exists near 
the edge at z < 30. Eventually, however, STb matches the background value at a large distance. The size of the region where 
we can detect STb becomes larger with increasing mass, and reaches almost 10 Mpc . 



4 THE ANGULAR POWER SPECTRUM OF 21 CM BRIGHTNESS TEMPERATURE DUE TO PBHS 

In order to study the potential of cosmological 21 cm observations to give a constraint on the PBH abundance, we evaluate 
the angular power spectrum of 21 cm fluctuations due to PBHs with mass M, 

C e (z obs ) = J dr W2( y obB) P(l/r,r), (17) 

where zi bs is the redshift by choosing the observation frequency v ob s = (1 + 2obs)f2i with V21 = 1420. MHz, r is the comoving 
radial distance, W(r, r b s ) is the window function of a 21 cm observation, and P(k, r) is the power spectrum of the brightness 
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Figure 4. The brightness temperature fluctuations as the function of the comoving radius. The left panel shows the dependence on 
the PBH's mass at z = 30. The right panel shows the dependence on the redshift for a PBH with M = 10 3 Mq. In both panels, the 
representations of lines are same as in Fig. [T] 



temperature fluctuations at r. In Eq. (|17|l . since we are only interested in the brightness temperature fluctuations on small 
scales (£ 10), we employ the Limber approximation. 

We obtain the power spect rum of the brightness temperature fluctuations P(k,r) by PBHs with mass M, following 
the halo formalism ( Seliak 2000). We assume that the spacial distribution of PBHs is totally random. Therefore, the power 
spectrum due to PBHs can be expressed only by the Poisson contribution, 

P(k) = n PBU {M)\ATi(k)\ 2 , (18) 

where upbh(M) is the comoving number density of PBHs with M. The number density of PBHs is described by using the 
density parameter flpBH as 

n PBH (M) = ^- = 1.36x10 (—jl^—j Mpc , (19) 



where p c is the critical density at the present epoch. In Eq. (1181) . ATb(k) is the Fourier transform of the brightness temperature 
fluctuations from the background value, 

AT b (k) = 4tt I x 2 dx(5T b (x) - 8T w ) 8m{ \ k \ (20) 
J xk 

where x is the comoving radius and STbo is the background brightness temperature given by Eqs. (£9|) and (|16[1 with the 
background kinetic temperature, Tko- 

Using 8Tb obtained in the previous section as shown in Fig. [4] we calculate the angular power spectrum. The results are 
shown in Fig. [5] Here we set J7pbh = 10 -11 . For simplicity, we assume that the window function has a strong peak at r — r b s 
and can be approximated as W 2 (r, r Q b s ) = 8(r — r D b s ). Since the promising signal of the cosmological 21 cm fluctuations is one 
due to t he primordial density flu ctuations, we also plot the angular power spectrum of those fluctuations obtained through 
CAMB jLewis fc Challinorll2007h . 

In the left panel of Fig. O it is shown that the peak amplitude of the angular spectrum is independent on PBH's mass. 
The peak locations shift toward large scales for massive PBHs. Moreover, the power law index of the spectrum matches the 
one due to the primordial density fluctuations on small £'s. The overall amplitude of the spectrum is simply proportional to 
f2pp.H or the PBH number density. For Qpbh = 10~ 12 , it is shown that the angular spectrum of PBHs with mass 
matches with the one due to the primordial density fluctuations at z = 30. It is clear that the spectrum with heavier (lighter) 
PBHs becomes comparable with primordial density fluctuations' one if Opbh is smaller (larger). To match these two spectra, 
we find the critical value of fipBH as ^pbhc = 10 _12 (M/10 3 Mq)~ ' 2 for z = 30. The PBH spectrum dominates over the one 
due to the primordial density fluctuations if ^pbh exceeds f^pBHc- 

The right panel of Fig.[S]shows that, although the peak location does not depend on the redshift, the redshift dependence of 
the peak amplitude is different from that due to the primordial density fluctuations. As the result, while the PBH contribution 
for fipBHc is subdominant at z = 40, even the PBHs with roughly 0.1 x JIpbhc can produce the spectrum comparable with 
that due to the primordial density fluctuations at z = 20. 

The brightness temperature profile at lower redshifts, z < 30, has both positive and negative peaks as shown in Fig. m 
The reason to have a negative peak is because there exists the region where the brightness temperature is lower than the 
background, which is rather difficult to see in the figure. Accordingly, the angular power spectrum due to PBHs also has two 
peaks at z < 30 as shown in Fig. [5] The peak on a larger scale is due to the negative peak of ATb, while that on a small scale 
is due to the positive one. 
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Figure 5. The angular power spectrum of 21 cm brightness temperature fluctuations produced by PBHs. comoving radius. The left 
panel shows the dependence on the PBH mass at z = 30. The dotted, solid and dashed lines represent for PBHs with M = 10 5 Mq, 
10 3 Mq and 10 Mq, respectively. The right panel shows the dependence on the redshift for PBHs with M = 10 3 Mq. The dotted, solid 
and dashed lines are for PBHs at z = 40, z = 30 and z = 20, respectively. We assume fippH = 10 -12 in both panels. For comparison, 
we also plot the angular power spectrum due to the primordial density fluctuations as thin lines. 



At the last of this section, we discuss the constraint by a future observation such as the square k ilometer arr ay (SKA)-like 
interferometer. The noise power spectrum of an observation including the beam effects is given by (|Knoxlll995l l 



£(£ + 1) 
2tt 



1(1 + 1) 


( D\T S y S ^ 




' 1(1 + 1)' 


iobsAi/ 




) exp 


P 



(21) 



where ^4 e fi/r s y S is the sensitivity (an effective area divided by the system temperature), t Q ^ s is the observation time, Av is 
the frequency bandwidth, D is the length of the baseline and lb is given by £t = 4Vln2/f?/ u , with the resolution 6f w ~ X/D. 



Here we adopt the current design of SKA0; 
design, the noise power spectrum is evaluated as 

i(t + l 



2000 m 2 K" 



^obs 



1000 hour, Av = 1 MHz and D = 5 km. In this 



2tt 



-N, 



73 mK' 



1000 



30 



exp 



z + 1 



2500 30 



(22) 



The noise power spectrum exponentially grows on larger multipoles than £i — 2500(30/(1 + z)). The amplitude of the angular 
power spectrum due to the PBHs is scaled by £2pbh. Accordingly we find that, in order to overdominate the noise spectrum, 
fipBH ~ 10" 9 (Af/10 3 Mq)-°- 2 is required for M > 1OM at z = 30 and Q PB h ~ 10" 10 (M/10 3 Mq)~ ' 2 is for M > 10M Q at 
z = 20. Since the angular spectrum due to PBHs with mass less than 10 Mq has a peak on larger multipoles than tb in the 
SKA design, there is no opportunity to measure the anisotropy spectrum due to such small mass PBHs by SKA. 



5 CONCLUSION 

We have investigated the 21 cm signal produced by massive PBHs whose masses are larger than 10 Mq. Assuming the power- 
law spectrum of X-ray photons from an accretion disk, we have studied the ionization and heating of IGM gas near a PBH 
and evaluated the differential 21 cm brightness temperature. We have shown that a PBH can induce an observable signal of 
differential 21 cm brightness temperature. The size of the region where we can find the differential brightness temperature 
typically reaches 1-10 Mpc for our interested PBH mass range while the size depends on the PBH mass. 

We have also calculated the angular power spectrum of 21 cm fluctuations due to PBHs. The peak position of the 
angular spectrum depends on the PBH mass, while the amplitude is independent. Comparing it with the angular power 
spectrum caused by primordial density fluctuations, we have found that both of them become comparable if JIpbh = 
10" 12 (M/10 3 M©) -0 - 2 at z = 30 and 10" 13 (Af/10 3 Mq)~ ' 2 at z = 20 for PBH's mass from 10 Mq to 10 8 Mq. If the 
density parameter is larger than these values, the angular power spectrum due to PBHs exceeds the one from primordial 
fluctuations and can be measured. In other words, we cannot set constraints on the PBH density parameter below these 
values from 21 cm observations. If we consider the sensitivity of the SKA-like observation, for example, we can detect the 
signal of PBHs up to S^pbh = 10~ 9 (Af/10 3 M )~°- 2 at z = 30 and 10~ 10 (M/10 3 Mq) -0 ' 2 at z = 20 for PBHs with mass 
from 10 2 Mq to 10 s Mq. 

The ionization of IGM due to PBHs with such density parameters does not affect the global reionization history of 

1 http://www.skatelescope.org/ 
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the universe since reionization from each PBH only covers a tiny patch of the universe. Unlike reionization from first stars, 
therefore, such reionization has little impact on C MB temperature an isotropies. Accordingly the PBH density parameter 
constrained from WMAP data, that is Opbh < 10 -7 ( Ricotti et al. 20081 ). is several order of magnitude larger than the value 
we obtained above. In other words, we can conclude that 21 cm fluctuation observations have a potential to probe the PBH 
abundance which is impossible to access by CMB observations. 

The most theoretical uncertainty in this model is the flux of photons due to the accretion to PBHs. In this paper, 
we assume that the X-ray photon flux amplitude is a tenth of the Eddington l uminosity for simplicity. The luminosity is 
considered to depend on the matter accretion rate to a PBH. Ricotti et al. ( 20081 ) have studied the luminosity for the Bondi- 
Hoyle accretion in details. Although the luminosity depends on the PBH mass and feedback effect on the ionization and 
temperature, they have shown that the luminosity for a PBH with M = 10 3 M e is roughly a hundredth of the Eddington 
luminosity at z > 20. In our model, the brightness temperature profile near a PBH at a certain redshift depends on only 
PBH's mass and the amplitude of the angular power spectrum is scaled by JIpbh- Therefore, if we assume that the amplitude 
of the X-ray photon flux is a hundredth of the Eddington luminosity, JIpbhc, which is the required density parameter for 
PBHs to dominate the 21 cm fluctuations due to primordial density fluctuations, is 10 _11 (M/10 4 Mq) -0 ' 2 at z — 30 and 
10" 12 (Af/10 4 Mo) -0 ' 2 at z = 20 for PBH's mass from 10 Mq to 10 s Mq. 
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